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Abstract.  Since  the  last  presentation  of  SC-3  on  numerical  values  of 
fundamental  geodetic  parameters  at  the  IAU  General  Assembly  at  Kyoto 
in  1997  there  were  some  conceptual  as  well  as  fundamental  numerical 
changes.  The  four  basic  parameters  of  geodetic  (ellipsoidal)  reference 
systems  (GRS)  can  no  longer  be  considered  as  constant  with  time:  J2, 
a ,  u,  and  GM  have  to  be  replaced  by  clearly  (±10-8  or  better)  specified 
mean  values  or  have  to  be  associated  with  a  specific  epoch  or,  in  case 
of  GM,  with  specific  reference  frames  (a  =  semi-major  axis  of  Earth 
ellipsoid,  J2  =  second  degree  zonal  harmonic  of  geopotential,  uj=  spin  of 
Earth  rotation).  In  case  of  (a,  J2  ....)  associated  tidal  reductions  must 
be  specifically  defined  in  view  of  particular  applications  and  significant 
differences  between  different  tidal  reduction  types.  Or  we  may  replace  “a” 
by  a  quantity  which  is  independent  of  tides  like  the  geopotential  at  the 
geoid,  Wo,  where,  however,  also  temporal  changes  are  now  discussed.  The 
official  geodetic  reference  systems  such  as  GRS  80  and  WGS  84  (revised  in 
97-form)  are  also  no  longer  truly  representing  reality;  a  new  system  GRS 
2000  is  desired.  We  are,  meanwhile,  able  to  define  and  determine  tidal  and 
non-tidal  (secular,  periodic,  aperiodic)  variations  of  some  fundamental 
geodetic  parameters.  Others  are  under  investigation.  New  precession 
and/or  nutation  formulas  to  be  adopted  by  IAU  in  2000  or  later  would 
imply,  again,  changes  in  geodetic  parameters  such  as  H  =  hydrostatic 
flattening.  Those  and  related  other  consequences  are  considered. 


1.  Introduction 

In  comparison  to  the  report  of  SC  3  presented  at  the  IUGG-General  Assembly  at 
Birmingham  (July  1999)  which  contains  an  overview  on  all  relevant  fundamental 
“constants”  of  geodesy  there  is  relatively  little  change  of  numerical  “current  best 
estimates”  so  that  the  list  of  updates  is  relatively  short. 

On  the  other  hand,  there  are  some  conceptual  changes  and  revisions  of  ac¬ 
curacies  of  such  “best  estimates”  which  have  to  be  taken  into  account.  However, 
in  agreement  jointly  with  IAU  those  official  parameters,  which  are  no  longer  up- 
to-date  and  which  in  the  past  were  sometimes  considered  as  constants,  as  well 
as  their  temporal  changes  have  to  be  considered  and  revised. 

Moreover,  remeasurements  of  the  “Newtonian  Gravitational  Constant”  G 
led  to  similar  conclusions  which  were  encountered  a  few  years  ago  with  the  spin 
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rate  of  the  Earth’s  rotation  uj.  In  case  of  u  we  had  to  truncate  the  official  value 
because  it  did  not  reflect  the  temporal  changes  sufficiently. 

In  case  of  G  we  are  now  more  sceptical  than  before  of  the  accuracy  of 
previous  determinations  of  “big  G ”  due  to  new  determinations  which  deviate 
substantially  from  the  previous  determinations.  However,  more  interesting  to 
astronomy  and  to  the  IAU  is  the  “geocentric”  gravitational  constant  GM  where 
M  is  the  mass  of  the  Earth.  GM  is  determined  from  satellite  orbits  using 
Kepler’s  third  law  and  is  one  of  the  best  determined  constants  (or  parameters, 
if  we  include  possible  temporal  changes). 

Quantum  physics  is  now  of  stronger  impact  to  geodesy  than  in  the  past.  For 
high  precision  time  measurements,  the  “cold  fountain”  principle  is  now  leading 
to  better  measurements  of  intervals  than  hydrogen  masers.  Moreover,  results 
of  macroscopic  absolute  gravimetry,  which  now  reach  accuracies  of  several  mi- 
crogals,  have  meanwhile  been  reconfirmed  by  “dropping  atoms”  methods  where 
atomic  interferomety  leads  to  similar  accuracies  under  laboratory  conditions. 
Consequently,  the  direct  measurement  of  “small  g”  i.e.,  gravity  at  the  surface  of 
the  Earth  as  vertical  derivative  of  the  geopotential  including  its  variations  with 
time  and  in  space,  will  in  the  future  be  more  reliably  determined.  Thus,  for  a 
conservative  held,  —~g=  VIE,  where  IE=geopotential.  Precise  time-keeping  is 
affected  by  this  fact  because  the  reduction  of  atomic  clock  measurements  to  the 
geoid  in  determing  TAI  involves  the  potential  difference  between  the  level  of  the 
observation  site  and  the  geoid.  The  exact  determination  of  the  geoid  height  of 
the  clock  location  is  no  longer  a  problem. 

It  is  important  to  realize  that  the  potential  IEo  at  the  geoid  is  indepen¬ 
dent  of  the  tidal  reductions  whereas  the  semi-major  axis  of  the  Earth’s  ellipsoid 
(in  the  best  fitting  sense)  depends  significantly  on  the  underlying  tidal  regime. 
As  different  tidal  systems  are  necessary  for  different  applications  and  purposes 
(which  is  not  always  well-understood  in  the  scientific  community)  we  are  con¬ 
fronted  with  various  (basically  three)  numerical  values  of  the  semi-major  axis 
of  the  Earth  ellipsoid  which  may  lead  to  confusion.  We  may  distinguish  “mean 
tide”,  “tide  free”  and  “zero  frequency  tide”  values. 


2.  A  new  Geodetic  Reference  System 

Meanwhile  global  satellite  altimetry  in  combination  with  continental  height  sys¬ 
tems  allows  us  to  determine  IEo  with  accuracy  comparable  to  the  accuracy  of 
the  semi-major  axis  a.  In  1980  the  consistent  ellipsoidal  (in  the  sense  of  Pizzetti- 
Somigliana)  Geodetic  Reference  System  was  derived  from  the  four  parameters 

(a,  GM,  J2,  u)  (1) 

where  J2  takes  the  role  of  the  flattening  /.  It  represents  the  second  zonal  har¬ 
monic  of  the  geopotential  and  is  basically  a  function  of  /;  moreover,  it  is  directly 
derived  from  satellite  orbit  analysis.  Theoretically,  these  four  parameters  are  in¬ 
dependent  of  each  other.  They  consistently  and  uniquely  determine  the  GRS  80 
including  its  external  spheroidal  gravity  field. 
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Bursa  et  al.  (1999)  have  meanwhile  determined  a  new  rounded  numerical 
value 

W0  =  (62636856.0  ±  0.5)m2S-2  (2) 

for  the  time  1993-1999.  No  significant  secular  trend  was  detected  for  the  mean 
sea  surface  (a  major  part  of  the  Earth’s  surface)  during  that  time.  All  major 
continental  height  systems  were  included  in  the  evaluation  besides  the  sea  surface 
using  Topex/Poseidon  data.  This  value  is  indeed  independent  of  the  underlying 
tidal  regime  whereas,  depending  on  the  specific  tidal  correction,  the  semi-major 
axis  varies  by  20  cm.  Depending  on  the  “inverted  barometer”  correction  to  be 
applied  to  the  sea  surface,  Wo  varies  by  0.33  m2s~2.  The  authors,  however, 
finally  preferred  to  omit  the  “inverted  barometer”  hypothesis. 

Grafarend  and  Ardalan  (1999)  computed  a  new  Geodetic  Reference  Datum 
GRD  2000  in  the  Somigliana-Pizzetti-sense  as  described  above.  However,  they 
adopted 

W0  =  (62636855.72  ±  0.5)m2s“2  (3) 

based  on  the  Finnish  height  datum  which  differs  by  only  1.8  cm  from  the  global 
geoid  (Bursa  et  al.,  1999).  If  we  assume  that,  due  to  systematic  perturbations, 
any  (absolute)  accuracy  of  better  than  5  cm  is  purely  formal  (Ries,  priv.  comm., 
2000)  in  satellite  altimetry  of  the  Topex-Poseidon  type,  a  deviation  of  1.8  cm 
can  be  ignored. 

It  appears,  consequently,  optimal  to  use  the  following  parameters  for  a  new 


GRS  2000: 

Wo  =  (62636  856.0T0.5)  m2s~ 2  (rounded)  (4) 

GM  =  (398  600  441.8  ±  0.8)  m3s~ 2  (Bursa  et  al,  1999)  (5) 

uj  =  7  292  115x  10-11  rad  s-1  (nominal)  (6) 

J2  =  (—4. 84165371736 X  10_4)±(3. 56  xUT11  )( “tide-free”),  EGM  96  (7) 


The  IAG  Executive  Committee  at  Birmingham  in  1999  considered  it  “pre¬ 
mature”  to  introduce  GRS  2000  now  but  in  view  of  the  new  system  to  be  adopted 
by  IAU  in  2000  or  2001  a  new  GRS  2000,  perhaps  based  on  updated  numerical 
values,  may  be  appropriate.  The  semi-major  and  semi-minor  axes  of  such  a  GRS 
2000  would  deviate  from  those  of  GRS  80  by  about  0.4  m. 

The  interrelation  between  these  four  parameters  is  relatively  weak;  Wq  de¬ 
pends  slightly  on  GM .  Other  parameters  such  as  the  hydrostatic  flattening  of  the 
Earth  can  only  be  well-determined  and  updated  when  new  nutation-precession 
parameters  have  been  adopted  by  the  IAU. 


3.  Conclusions 

For  astronomical  purposes  ellipsoidal  models  appear  sufficient  so  the  discussion 
of  higher- degree  approximations,  as  presently  considered  in  geodesy,  is  not  nec¬ 
essary  here. 

A  fundamental  problem  lies  in  the  fact  that  GM ,  a  and  Wo  are  more  or 
less  interrelated  (pseudo-observations)  and  not,  as  assumed  in  theory  of  the 
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Somigliana-Pizzetti  concept,  strictly  independent.  When  Wo  was  introduced 
as  a  primary  constant,  an  improved  determination  (by  satellite  altimetry,  tide 
gauge  networks  and  repeat  GPS)  would  be  needed  in  order  to  verify  the  present 
(1999)  best  estimate  of  Wo=  (162636856. 3±0. 5)  m2s~ 2  primarily  based  on  the 
determination  of  J.  Ries.  The  zero-tide  and  tide-free  solutions  appear  to  be  the 
most  appropriate  for  astrometry,  in  general. 

Consequently,  our  latest  list  of  fundamental  parameters,  reconciled  with 
IERS  Conventions  in  1999,  now  needs  only  little  and  minor  revisions  of  the 
numerical  values,  but  rather  numerous  specifications. 

The  official  introduction  of  a  Geodetic  Reference  System  “GRS  2000”  is 
strongly  recommended  in  view  of  the  significant  deviations  of  up-to-date  ellip¬ 
soidal  gravity  fields  of  the  Earth  from  GRS  80  as  an  official  reference  system. 
An  obstacle  to  GRS  2000  is  the  fact  that  only  recently  WGS  84  was  officially 
adopted  by  international  navigation  authorities  (air  and  sea  navigation)  where 
WGS  84  is  almost  identical  with  GRS  80.  The  updated  list  of  ‘constants”  pre¬ 
sented  at  Birmingham  in  1999  is  given  in  appendices  A  and  B. 
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Appendix  A 

A1  Current  (2000)  best  estimates  of  the  parameters  of  common  rel¬ 
evance  to  astronomy,  geodesy,  and  geodynamics 

SI  units  are  used  throughout  (with  the  exception  of  the  TDB-value  below  (4)). 
(Si-values  can  be  associated  with  TCB  or  TCG). 

-  velocity  of  light  in  vacuum 

c  =  299792458  ms-1  (1) 

-  Newtonian  gravitational  constant 

G  =  (6672.59  ±  0.30)  X  10-u m3s~2 kg-1  .  (2) 

-  Geocentric  gravitational  constant  (including  the  mass  of  the  Earth’s 
atmosphere);  reconfirmed  by  J.  Ries  (1998,  priv.  comm.) 


GM  =  (398600441.8  ±  0.8)  X  106m3s“2. 


(3) 
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For  the  new  EGM  96  global  gravity  model  GM  =  398600441.5  X  106m3s  2  was 
adopted. 

In  TT  units  (Terrestrial  Time)  the  value  is 

GM  =  (398600441.5  ±  0.8)  X  106m3s“2.  (4) 

Note  that  if  expressed  in  old  TDB  units  (solar  system  Barycentric  Dynamical 
Time),  the  value  is 

GM  =  398600435.6  X  106m3s“2. 

Based  on  well-known  transformation  formulas  we  may  relate  GM  in  Si-units 
to  TT/TCG/TCB;  see  IERS  Convention  1996  p.  85.  The  well-known  secular 
term  was  not  originally  included  in  the  GM(E)-analysis.  Therefore  it  was  related 
to  TT,  not  to  SI  or  (TCG,  TCB).  As  satellite  analyses  still  occur  without  the  sec¬ 
ular  term,  GM(E)  in  TT  is  still  of  geodetic  interest.  GM(E)  =  GM  of  the  Earth. 

-  Mean  angular  velocity  of  the  Earth’s  rotation 

w  =  7292115  X  lQ-^rads-1.  (5) 


Table  1.  Mean  angular  velocity  of  the  Earth’s  rotation  1978-1999. 


Year 

u)  [10  11rad  s  x] 

Year 

u)  [10  11rad  s  x] 

Mean  LOD  [ms /day] 

min:  1978 

7  292  114.903 

1995 

7  292  114.952 

- 

max:  1999 

292  115.063 

1996 

7  292  114.992 

- 

1997 

7  292  114.991 

- 

1998 

7  292  115.031 

1.37 

1999 

7  292  115.063 

0.99 

-  Long-term  variation  in  u 

duj/dt  =  (—4.5  ±0.1)  X  10-22rad  s~2 .  (6) 

This  observed  average  value  is  based  on  two  actual  components: 
a)  due  to  tidal  dissipation 

(d^/dt)tidal  =  (-6-1  ±  °-4)  x  10“22rad s“2.  (7) 

This  value  is  commensurate  with  a  tidal  deceleration  in  the  mean  motion  of 
the  Moon  n 


diif  dt  =  (  —  25.88  ±  0.5)  arcsec  cy  2 . 
b)  non-tidal  in  origin 


(8) 
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(d(jJ / dt)non-tidal  =  (  +  1'6  ±  °'4)  X  10  22 rad  S  2-  (9) 

-  Second-degree  zonal  geopotential  (Stokes)  parameter  (tide-free,  fully 
normalized,  Love  number  £;2  =  0.3  adopted),  in  agreement  with  EGM  96, 

72  =  4.84165371736  X  10“4  ±  3.56  X  10“n.  (10) 

To  be  consistent  with  the  I.A.G.  General  Assembly  Resolution  16,  1983 
(Hamburg),  the  indirect  tidal  effect  on  J2  should  be  included.  Then  in  the 
zero-frequency  tide  system  (JGM-3) 

J2  =  (1082635.9  ±  0.1)  X  10“9.  (11) 


Table  2.  The  Stokes  second-degree  zonal  parameter;  marked  with  a 
bar:  fully  normalized;  £;2  =  0.3  adopted  for  the  tide-free  system. 


Geopotential  model 

Zero-frequency  tide  system 

72[10“6]  «/2[10-6] 

Tide-free 

72[10“6]  «/2[10-6] 

JGM-3 

484.16951  1082.6359 

484.16537  1082.6267 

EGM  96 

484.16537 

-  Long-term  variation  in  J2 

dJ2/dt  =  -(2.6  ±  0.3)  X  10-9cr/-1.  (12) 

-  second-degree  sectorial  geopotential  (Stokes)  parameters  (conventional,  not 
normalized,  geopotential  model  JGM-3) 

Jf  =  (1574.5  ±  0.7)  X  10-9,  (13) 

Si  =  -  (903.9  ±  0.7)  X  10“9,  (14) 

J2, 2  =  ^[(./l)2  +  (Si)2]  =  (1815.5  ±  0.9)  X  10“9.  (15) 


Table  3.  The  Stokes  second-degree  sectorial  parameters. 

marked  with  a  bar:  fully  normalized 


Geopotential  model 

c2  [io-6] 

^2  [10-6] 

JGM-3 

EGM  96 

2.43926 

2.43914 

-1.40027 

-1.40017 
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Only  the  last  decimal  is  affected  by  the  standard  deviation. 

For  EGM  96  Marchenko  and  Abrikosov  (1999)  found  more  detailed  values: 


Table  4.  Parameters  of  the  linear  model  of  the  potential  of  2nd  degree. 


Harmonic  coefficient 

Value  of  coefficien  xlO  b 

Temporal  variation  xlO LL[yr  i] 

C  20  =  —  J  2 

-484.165371736 

1.16275534 

C21 

-0.00018698764 

-0.32 

S  21 

0.00119528012 

1.62 

c22  =  -4 

2.43914352398 

-0.494731439 

S  22 

-1.40016683654 

-0.203385232 

-  Coefficient  H  associated  with  the  precession  constant 

H  =  (C  —  0.5  (A  ±  B  j)  /C  =  (3273763  ±  20)  X  10“9.  (16) 


-  The  geoidal  potential  ITo  and  the  geopotential  scale  factor  Ro  =  GM/Wo 
recently  derived  by  Bursa  et  al.  (1998)  read 


W0=  (62  636  855.611T  0.5)  m2s~ 2, 
Ro  =  (6  363  672.58  ±0.05)  to, 


(17) 


Wq  =  (62636856.4  ±0.5)  m2s  2  J.  Ries  (priv.  comm,  1998)  found  globally. 


If  Wq  is  preserved  as  a  primary  constant  the  discussion  of  the  ellipsoidal 
parameters  could  become  obsolete  as  the  Earth  ellipsoid  is  basically  an 
artifact.  Modelling  of  the  altimeter  bias  and  various  other  error  influences 
affect  the  validity  of  ITo-determination.  The  variability  of  Wo  and  Ro  was 
studied  by  Bursa  et  al.  (1998)  recently;  they  detected  interannual  variations  of 
ITo  and  Rq  amounting  to  2  cm. 


The  relativistic  corrections  to  ITo  were  discussed  by  Kopejkin  (1991).  See  his 
formulas  (67)  and  (77)  where  tidal  corrections  were  included.  Whereas  he 
proposes  average  time  values,  Grafarend  insists  in  corrections  related  to 
specific  epochs  in  order  to  illustrate  the  time-dependence  of  such  parameters  as 
ITo,  Gill,  Jn,  which  are  usually,  in  view  of  present  accuracies,  still  treated  as 
constants  in  contemporary  literature. 


Based  on  recent  GPS  data,  Grafarend  and  Ardalan  (1997)  found  locally  (in  the 
Finnish  Datum  for  Fennoscandia)  ITo  =  (6  263  685. 58±  0.36)  kgalm. 

The  temporal  variations  were  discussed  by  Wang  and  Kakkuri  (1998),  in 
general  terms. 


Mean  equatorial  gravity  in  the  zero-frequency  tide  system 
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ge  =  (978032.78  ±  0.2)  x  10“5ms-2. 

(18) 

-  Equatorial  radius  of  the  Reference  Ellipsoid  (mean  equatorial  radius  of  the 
Earth)  in  the  zero-frequency  tide  system  (Bursa  et  al. ,  1998) 

a  =  (6378136.62  ±  0.10)  to. 

(19) 

-  The  corresponding  value  in  the  mean  tide  system  (the  zero- 
and  indirect  tidal  distortion  included)  comes  out  as 

frequency  direct 

a  =  (6378136.72  ±  0.10)  m  , 

(20) 

and  the  tide-free  value 

a  =  (6378136.59  ±  0.10)  m  . 

(21) 

The  tide  free- value  adopted  for  the  EGM-96  gravity  model  reads 
a  =  6378136.3  m. 

-  Polar  flattening  computed  in  the  zero-frequency  tide  system,  (adopted  GM, 


uj,  and  J 2  in  the  zero-frequency  tide  system) 

l//  =  298.25642  ±  0.00001.  (22) 

The  corresponding  value  in  the  mean  tide  system  comes  out  as 
l//  =  298.25231  ±  0.00001,  (23) 

and  the  tide-free 

l//  =  298.25765  ±  0.00001.  (24) 

-  Equatorial  flattening  (geopotential  model  .JGM  —  3) 

1/oq  =  91026  ±  10.  (25) 

-  Longitude  of  major  axis  of  the  equatorial  ellipse,  geopotential  model  JGM  —  3 

Aa  =  (14.9291°  ±  0.0010°)  W.  (26) 

In  view  of  the  small  changes  (see  Table  3)  of  the  second  degree  tesserals  it  is 


close  to  the  value  of  EGM  96.  We  may  raise  the  question  whether  we  should 
keep  the  reference  ellipsoid  in  terms  of  GRS  80  (or  an  alternative)  fixed  and 
focus  on  Wo  as  a  parameter  to  be  essentially  better  determined  by  satellite 
altimetry,  where  however  the  underlying  concept  (inverted  barometer, 
altimeter  bias  etc.)  has  to  be  clarified. 

-  Coefficient  in  potential  of  centrifugal  force 
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Table  5.  Equatorial  flattening  oq  and  Aa  of  major  axis  of  equatorial 
ellipse. 


Geopotential  model 

l/«l 

Aa  [deg.] 

JGM  -  3 

91026 

14.9291  W 

q  =  ( uj2ci3)/(GM )  =  (3461391  ±  2)  X  10“9.  (27) 

Computed  by  using  values  (3),  (5)  and  a  =  6378136.6. 

-  Principal  moments  of  inertia  (zero-frequency  tide  system),  computed  using 
values  (11),  (15), (3), (2)  and  (16) 


(C  -  A) /(Mai )  =  -h  +  2/2,2  =  (1086.267  ±  0.001)  X  10“6, 

(C  -  B)/(Ma2)  =  J2-  2/2,2  =  (1079.005  ±  0.001)  X  10“6, 

( B  -  A) /(Mai)  =  4/2,2  =  (7.262  ±  0.004)  x  10“6 

(28) 

Mai  =  ( GMa2)/G  =  (2.43014  ±  0.00005)  X  10 38kgm2, 

(a0  =  6378137  m)  ; 

(29) 

C  -  A  =  (2.6398  ±  0.0001)  X  1035  kg  m2 , 

C-B  =  (2.6221  ±  0.0001)  X  1035  kg  m2 , 

B  -  A  =  (1.765  ±  0.0001)  X  1033  kg  m 2  ; 

(30) 

C/(Mal)  =  J2/H  =  (330701  ±  2)  x  10“6  , 

A/ (Mai)  =  (329615  ±  2)  x  10“6  , 

(31) 

B/(Ma §)  =  (329622  ±  2)  x  10“6  ; 

A  =  (8.0101  ±  0.0002)  x  1037  kg  m2  , 

(32) 

B  =  (8.0103  ±  0.0002)  x  1037  kg  m 2  , 

A  =  (8.0365  ±  0.0002)  x  1037  kg  m2  , 

a  =  (C  —  B)/A  =  (327353  ±  6)  X  10“8  , 

7  =  (B  -  A)/C  =  (2196  ±  6)  X  10“8  , 

f3  =  (C  -  A)/B  =  (329549  ±  6)  X  10“8  . 

(33) 
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A2  Primary  geodetic  Parameters,  discussion 

It  should  be  noted  that  parameters  a,  /,  J2 ,  ge ,  depend  on  the  tidal  system 
adopted.  They  have  different  values  in  tide-free,  mean  or  zero-frequency  tidal 
systems.  However,  Wq  and/or  Ro  are  independent  of  tidal  system  (Bursa, 
1995).  The  following  relations  can  be  used: 

a(mean)  =  a(tide  —  free)  +  0.5(1  +  ks)  Ro  (SJ2)/ks  ,  (34) 

a(mean)  =  a(tide  —  free)  +  1.5(1  +  ks)  (SJ2)/ks  . 

a(zero  —  frequency)  =  a(tide  —  free)  +  0.5  R0SJ2,  (35) 

a(zero  —  frequency)  =  a(tide  —  free)  +  1.5  SJ2. 

ks  =  0.9383  is  the  secular  Love  number,  (SJ2  is  the  zero-frequency  tidal 
distortion  in  .J2).  First,  the  internal  consistency  of  parameters  a,  ITo,  (Ro)  and 
ge  should  be  examined: 

(i)  If 

a  =  6378136.7  m 

is  adopted  as  primary,  the  derived  values  are 
Wo  =  62636856.88  m2s~2  , 

(Ro  =  6363672.46  m)  , 
ge  =  978032.714  x  10“5  ms~2 

(ii)  If 

Wo  =  (62636855.8  ±  0.5)  m2s~2  , 

Ro  =  (6363672.6  ±  0.05)  m  , 

is  adopted  as  primary,  the  derived  values  are  (mean  system) 
a  =  6378136.62  m, 
ge  =  978032.705  X  10“5  ms~2 . 

(iii)  If  (18) 

ge  =  (978032.78  ±  0.2)  x  10“5  ms~2 , 


(36) 

(37) 

(38) 


(39) 

(40) 
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is  adopted  as  primary,  the  derived  values  are 

a  =  6378136.38  m,  (41) 

W0  =  62636858.8  m2s~ 2,  (42) 

Ro  =  6363672.26  m.  (43) 


There  are  no  significant  discrepancies,  the  differences  are  about  the  standard 
errors. 

However,  the  inaccuracy  in  (iii)  is  much  higher  than  in  (i)  and/or  (ii).  That  is 


why  solution  (iii)  is  irrelevant  at  present. 

If  the  rounded  value 

W0  =  (62636856.0  ±  0.5)  m2s~ 2  ,  (44) 

Ro  =  (6363672.6  ±  0.1)  m  (45) 

is  adopted  as  primary,  then  the  derived  length  of  the  semi-major  axis  in  the 
mean  tide  system  comes  out  as 

a  =  (6378136.7  ±  0.1)  to,  (for  zero-tide:  6378136.6),  (46) 

which  is  just  the  rounded  value  (20),  and  (in  the  zero  frequency  tide  system) 

=  (978032.7 ±  0.1)  X  10“5tos-2.  (47) 


However,  SC  3  recommends  that,  at  present,  GRS  1980  should  be  retained  as 
the  standard. 


A3  Consistent  set  of  fundamental  constants  (1997) 

-  Geocentric  gravitational  constant  (including  the  mass  of  the  Earth’s 
atmosphere) 

GM  =  398600441.8  ±  0.8)  X  106to3s“2  , 

[value  (3)] 


-  Mean  angular  velocity  of  the  Earth’s  rotation 

w  =  7292115  X  10“n  rads-1 
[value  (5)] 
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-  Second-degree  zonal  geopotential  (Stokes)  parameter  (in  the  zero-frequency 
tide  system,  Epoch  1994) 

■J2  =  (1082635.9  ±  0.1)  X  10“9 
[value  (11)] 

-  Geoidal  potential 

W0  =  (62636856.0  ±  0.5)  m2s~2  , 

[value  (44)] 

-  Geopotential  scale  factor 

R0  =  GM /Wq  =  (6363672.6  ±  0.05)  m 
[value  (45)] 


-  Mean  equatorial  radius  (mean  tide  system) 

a  =  (6378136. 7±  0.1)  m 
[value  (46)] 

-  Mean  polar  flattening  (mean  tide  system) 

1  //  =  298.25231  ±  0.00001 
[value  (23)] 

-  Mean  equatorial  gravity 

ge  =  (978032.78  ±  0.1)  X  10“5  ms'2, 

[value  (18)]. 

Grafarend  and  Ardalan  (1999)  have  evaluated  a  (consistent)  normal  held  based 
on  a  unique  set  of  current  best  values  of  four  parameters  (W0,  oj,  -J2  and  GM) 
as  a  preliminary  “follow-up”  to  the  Geodetic  Reference  System  GRS  80.  It  can 
lead  to  a  level-ellipsoidal  normal  gravity  held  with  a  spheroidal  external  held  in 
the  Somigliana-Pizetti  sense.  By  comparing  the  consequent  values  for  the 
semi-major  and  semi-minor  axes  of  the  related  equipotential  ellipsoid  with  the 
corresponding  GRS-80  axes  (based  on  the  same  theory)  the  authors  end  up 
with  axes  which  deviate  by  -40  and  -45  cm,  respectively  from  GRS  80  axes  and 
within  standard  deviations  from  the  current  values  such  as  in  (21);  but  no 
g-values  are  given  until  now. 


Appendix  B 

B1  Zero-frequency  tidal  distortion  in  -J2 

(J2  =  —C2o) 


Report  oflAG  Special  Commission  3 


349 


SJ2  =  ks[(GMl)/(GM)\  (R/AeL)3  (. ft/ao )2  (E2  +  S2L)  + 

Ks[(GMs)I(GM)}  (r/Aqs)3  (R/oo)2  (E2+S2s)  , 

E2  =  —0.5  +  0.75  sin2  e0  , 

hn  =  |  (sin2  iL  -  ej)  +  \e\  (sin2  e0  -  sin2  iL)  , 
hs  =  -|e|  (l  -  |  sin2  e0)  , 

R  =  Ro  (ff^g  -  t772^2)  / 

GMl  =  4902.799  X  109m3s-2  (selenocentric  grav.  Const.), 

GMS  =  13271244.0  X  1013m3s“2, 

A®/,  =  384400  km  (mean  geocentric  distance  to  the  Moon), 

A®5  =  1  AU  =  1.4959787  X  1011  m  , 

a0  =  6378137  m  (scaling  parameter  associated  with  J2), 
e0  =  23°26,21"4  (obliquity  of  the  ecliptic), 

=  0.05490  (eccentricity  of  the  orbit  of  the  Moon), 
iL  =  5°0l9  (inclination  of  Moon’s  orbit  to  the  ecliptic), 
eg  =  0.01671  (eccentricity  of  the  heliocentric  orbit  of  the  Earth-Moon 
bary  center), 

v  =  a0/R0  =  1.0022729; 

ks  =  0.9383  (secular-fluid  Love  number  associated  with  the  zero-frequency 
second  zonal  tidal  term); 

&  J2  =  —SC2o  =  (3.07531  X  10-8)  ks  (conventional); 

& J 2  =  —SC 20  =  (1.37532  X  10-8)ks  (fully  normalized). 

L  =  lunar 
S  =  solar 


B2  Definition 

Because  of  tidal  effects  on  various  quantities,  the  tide-free,  zero-frequency  and 
mean  values  should  be  distinguished  as  follows: 

-  A  tide-free  value  is  the  quantity  from  which  all  tidal  effects  have  been 
removed. 

-  A  zero-frequency  value  includes  the  indirect  tidal  distortion,  but  not  the 
direct  distortion. 

-  A  mean  tide  value  includes  both  direct  and  indirect  permanent  tidal 
distortions. 

Acknowledgments.  This  report  is  basically  an  updated  version  of  M. 
Bursa’s  SC  3  report  presented  in  1995  with  some  new  material  added. 
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